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Abstract 

I discuss the constraints on the lattice spacing, a, the quark masses, to, the box size, L, and 
particularly the residual mass, m res , such that one can successfully calculate phenomenologically 
interesting quantities using Domain Wall fermions (DWF). The constraints on a, m, and L are 
largely common with other improved fermion discretizations, and I emphasize that the improved 
chiral symmetry of DWF does not remove the need for simulations with a significant range of 
lattice parameters. Concerning m res , I compare the analysis of chiral symmetry breaking to that 
with Wilson fermions, emphasizing that DWF are better than simply Wilson fermions with each 
chiral symmetry breaking effect reduced by a common factor. I then discuss the impact of non-zero 
m res both on generic hadronic quantities, and on matrix elements which involve mixing with lower 
dimension operators. 
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I. INTRODUCTION 



Simulations with domain- wall fermions (DWF) |2| have made impressive progress over 
the last few years. Algorithmic improvements (particularly the work of Ref. (3|) have allowed 
simulations with 2 + 1 flavors of dynamical fermions to begin in earnest. Precision results 
for important quantities have been obtained, e.g. Bx [H and the K —>■ n form factor 
albeit at a single lattice spacing. Thus it is a good time to look ahead a few years and study 
the strengths and limitations of the approach [6j]. 

The key attraction of DWF compared to Wilson-like fermions is the reduction in the size 
of chiral symmetry breaking. This comes at the cost of adding an extra flavor-like dimension, 
with concomitant increase in computational cost. The length of the extra dimension, L5, is 
chosen as a compromise between speeding up the simulations (by reducing L 5 ) and decreasing 
the size of chiral symmetry breaking (by increasing L 5 ). The main focus of this talk is 
to study the impact of the residual chiral symmetry breaking (XSB) on calculations of 
phenomenologically important quantities. 

I cannot resist setting things up as a choice between three alternatives: 

• DARN GOOD: No practical barrier to calculating some quantities of interest (spec- 
trum, Bk • • • ) with desired precision in next 5 years. 

• WONDERFUL: No practical barrier to calculating many quantities of interest (spec- 
trum, Bk, some K — > tttt matrix elements, . . . ) with desired precision in next 5 
years. 

• FANTASTIC: No practical barrier to calculating most quantities of interest (spec- 
trum, Bk, e' /e, condensate . . . ) with desired precision in next 5 years. 

There are no bad (or ugly 0]) options here — i.e. it is expected that the theory has the cor- 
rect continuum limit. The only theoretical concern is whether the effective four-dimensional 
action is local, and this can be checked numerically, as has been done for present simula- 
tions jl] . There is, however, the practical issue of how far calculations can be pushed given 
available computational resources. What is needed is an ensemble of lattice ensembles with 
a range of lattice parameters (a, rri£, m s , and box size L) comparable to that of the present 
MILC staggered-fermion ensemble. This would allow validation of the methods by compar- 
ing to experimentally measured quantities, independent predictions of those quantities which 
have been accurately calculated with staggered fermions, 1 and first controlled calculations 
of the quantities which are difficult to obtain using staggered or Wilson-like fermions. It is 
the latter quantities that are my focus here. 

The outline of this talk is as follows. In the next section I briefly recall the range of 
simulation parameters, in particular the light quark masses, that will likely be needed to 
extract precision results. This is to emphasize the importance of creating the ensemble of 
ensembles mentioned above, and to bring out a couple of points that are sometimes forgotten. 
I then turn to the core of this talk, namely an analysis of the impact of the residual XSB 
on various important quantities: pion properties, the quark condensate, electroweak matrix 
elements with and without power divergences, and the pion electromagnetic mass splitting. I 



Here I am injecting my opinion that it is plausible that rooted-staggered fermions give rise to the correct 
continuum limit, based on the considerations and work of others summarized in Refs. 
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also discuss the extent to which one can think of DWF as "Wilson-fermions-lite, " i.e. Wilson 
fermions with the XSB reduced (very significantly) in magnitude. I close by attempting a 
choice from the three options given above. 

I should stress that I am mainly pulling together work in the literature and attempting to 
provide a uniform discussion, although I am adding some new observations along the way. 

II. HOW SMALL DO a, m AND l/L NEED TO BE? 

In this section I make a few remarks on the range of simulation parameters that one 
should aim for. 2 This range obviously depends on the quantity under consideration and 
the desired precision. What I have in mind here are "basic" quantities (spectrum, decay 
constants, matrix elements without subtractions, . . . ) and percent-level accuracy. Nothing 
I say here is special to DWF, but conversely one needs to keep in mind that the positive 
features of DWF do not exempt them from these standard requirements. 

A. How small does a need to be? 

Although residual XSB introduces discretization errors linear in the lattice spacing a (at 
least for generic quantities, as discussed further below), the dominant such error for most 
quantities is of relative size 0(aA) 2 . In other words, DWF are (almost) automatically 0(a) 
improved. This puts DWF simulations in essentially the same situation as almost all the 
other large-scale simulations, i.e. those using overlap fermions, non-perturbatively 0(a) 
improved Wilson fermions, twisted-mass fermions, or staggered fermions. 

The scale characterizing discretization errors should be comparable to Aqcd ~ 0.3 GeV, 
and I start by taking A = 0.1 — 0.5 GeV. With this, the relative size of the leading dis- 
cretization errors is 



Since one has, or aims to have, two or more lattice spacings, this leading error should be 
largely removed by extrapolation to a = 0, leaving a significantly smaller residue. Thus it 
is not unreasonable to think that percent-level accuracy should be attainable with a min « 
0.1 fm. This is consistent with the RBC/UKQCD plans described at this meeting, which 
involve moving from the present runs with a « 0.12 fm to runs (already underway) at 
a « 0.09 fm. 

Is this assessment correct? Are simulations with a m i n 0.1 fm sufficient? I raise this 
issue because there is evidence that smaller lattice spacings are needed for improved Wilson 
fermions, despite the fact that the leading discretization errors are of the same parametric 
size as for DWF. 3 Some of the evidence is collected in Sommer's recent lectures He 
quotes the following examples, all using non-perturbatively improved Wilson fermions: 



2 I do not discuss constraints on L5 in this section — these will be discussed at length below. For the 
discussion in this section, and in particular for the "basic" quantities under consideration, I assume that 
L5 is fixed at a large enough value that residual XSB effects are smaller than other errors. 
Improved staggered fermions are also pushing to smaller lattice spacings — a m 0.06 fm at present, and 
possibly a w 0.045 fm in the future. The need for smaller lattice spacings is greater for staggered fermions, 
however, because of the need to disentangle taste breaking from the chiral limit. 




(1) 
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10-15% discretization errors at a = 0.1 fm seen in quenched m q and fx [ll|- These 
errors are larger than the range above, perhaps indicating that one should use scales 
as large as A = 750 MeV. 

Similar effects seen for Nf = 2 12l . 

Large differences between results using different definitions of Za (which should agree 
in the continuum limit) for a > 0.1 fm. 13 . 



Perhaps because of this, recent high-precision calculations with Wilson-like fermions have 
worked with lattice spacings down to a min w 0.05 fm [14], |l5 . 

One should also note, however, that improved Wilson fermions do differ from DWF. In 
particular, the leading sub-dominant error for improved Wilson fermions is proportional to 
a 3 , an error which is much suppressed for DWF because it is proportional to the residual 
%SB. It could be, therefore, that discretization errors are smaller with DWF, as has been 



observed for some quantities using overlap fermions [16|, |17||. Nevertheless, it seems prudent 
to expect larger discretization errors than the estimate in ([1]) suggests, at least for some 
quantities, and thus that one may need to work at significantly smaller lattice spacings than 
0.1 fm. 



B. How small does m q need to be? 



Concerning the light quark masses, it probably bears repeating that extrapolations to 
the physical value of the average light quark mass, mg = (m u + m^)/2 ps m^ hys /27, require 
that one work at small enough to be able to see and fit the non-analytic contributions. 
Examples of such "chiral logarithms" (taken from Ref . 18] , where further discussion may be 
found) are shown in Fig. [TJ Here I have plotted the results of continuum XPT taking mf ys = 
0.08 GeV, / ps 0.09 GeV, and, for the Gasser-Leutwyler coefficients, L 5 = 1.45 x 10~ 3 , 
L$ = 10~ 3 and L 4 = L 6 = at a scale fi = m p [l9|. These curves are indicative of what one 
should find when doing unquenched chiral extrapolations (and not precise predictions). 




FIG. 1: Examples of curvature due to the presence of chiral logarithms. m s is held fixed, while mg 
is varied down to its physical value. See text for more discussion. 

My conclusion is that one needs to have accurate results down to m^/mP hys ps 0.1 for 
precision extrapolations. This is of course a generic conclusion — individual quantities or 
ratios may have smaller non-analytic terms and be more easily extrapolated [20j. It is also 
"old news" — the MILC collaboration has worked down to such masses and found that they 
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do allow precision extrapolations [2l|], and other simulations are aiming at, or have already 
attained, such small masses. But perhaps it is good to keep being reminded of it! 



C. How large does L need to be? 

Given the cost of simulating with DWF, work to date has been on relatively small lattices 
with a physical size of about 2 fm. Here I want to emphasize the importance of using larger 
lattices, given that finite volume effects can be substantial. I illustrate this in Fig. [2J which 
shows that at irti/m s = 0.1, finite volume shifts are at the few percent level. This is well- 
known, and indeed generation of L = 3 fm DWF lattices at a « 0.12 fm is underway in order 
to study, and reduce, finite volume effects. 
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FIG. 2: Example of finite volume effects, as predicted by one-loop xPT, assuming a = 0.1 fm, 
and taking L = 2.4fm (thick line-at bottom), 3.2fm (medium line-middle) and oo (thin line-top). 
Parameters as for other figures, except that / = 0.08 GeV. See text for more discussion. 

What is perhaps less well-known is that one-loop XPT does not, in general, give an 
accurate estimate of volume corrections — one needs to include (at least the dominant part 
of) the two-loop contributions [22]. To say it differently, one-loop results give, in general, a 



good indication of the magnitude of the effect but are not accurate enough to do a correction 
accurate at the percent level. Since two-loop calculations are hard, and are likely to be 
available for only a few quantities, a conservative conclusion is that one must work in volumes 
large enough so that the one-loop correction is smaller than the desired accuracy. This pushes 
one to even larger volumes, although how large depends on the quantity and the desired 
precision. 



III. IMPLICATIONS OF RESIDUAL CHIRAL SYMMETRY BREAKING 

Domain-wall fermions work by having left- and right-handed quarks (with small four- 
dimensional momentum) bound to different defects, in this case the two walls at either end of 
the fifth dimension. Chiral symmetry — the separate rotation of left- and right-handed quark 
fields — is broken by the overlap of the bound-state "wavef unctions". This rough description 
can be made more precise using a transfer-matrix formalism to describe propagation in the 
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fifth dimension [23], 24]. The result is that, for the interactions of the low-momentum four- 
dimensional quark fields (i.e. those with p<l/o), the size of XSB can be parameterized by 
the residual mass (the various definitions of will be discussed below) whose dependence on 



L 5 is schematically |8j, 125 



-OlLk 



+ 



P(0) 



(2) 



L5 L$ 

I take m res to be dimensionless in this talk, primarily because it simplifies the comparison 
with Wilson-like fermions. Some intuition concerning the form of m res will be useful in the 
following, so I give a brief discussion. 

The first term on the r.h.s. corresponds to the picture of modes exponentially localized on 
the walls, which thus have an overlap falling exponentially with L5. In the transfer matrix 
language, this term arises from delocalized eigenmodes, which occur for eigenvalues above 
the "mobility-edge" a. In perturbation theory a ~ 0(1), and depends on the negative mass 
used in the fifth dimension and on the choice of action [26l . |27 . 

The power-law term in ([2]) is not visible in perturbation theory, but arises in the transfer- 
matrix analysis from near-zero modes of the corresponding Hamiltonian, of which p(0) is the 
density per unit four- volume. A key observation is that the corresponding (four-dimensional) 
eigenf unctions are localized at the scale of the lattice spacing [28]. Crudely speaking, a 
left-handed quark propagating on one wall has, generically, only an exponentially damped 
amplitude to "tunnel" to the other wall (and thus become right-handed), but occasionally 
encounters a near-zero mode, which can transport it to the other wall with an essentially 
unsuppressed amplitude. There are, however, constraints on simultaneous "transports" of 



multiple fermions which will be important later |29J. As for a, the zero- mode density p(0) 
can be altered, and in particular reduced, by tuning the gauge and DWF actions. 

It is important to keep in mind that both delocalized and localized eigenmodes of the 
transfer matrix lead to a local effect in the long-distance four-dimensional effective theory, 
and in particular to a local mass term at leading order in the Symanzik expansion. Here 
local means exponentially localized at the lattice scale. For the localized eigenmodes, this 
is automatic. For the delocalized eigenmodes, however, this must arise from the collective 
effect of the modes in the vicinity of the mobility edge. Heuristically, the argument for 
locality is that, if the mobility edge is of 0(1/ a), i.e. of order its perturbative value, then, 
in the viewpoint in which the fifth dimension is a flavor index, one has L 5 four-dimensional 
flavors coupled with a mass matrix with generic entries of 0(1/ a). Diagonalizing this matrix 
leads to the single light Dirac fermion and L5 — 1 heavy fermions with masses of 0(1/ a). 
Integrating the latter out will lead to a local four dimensional theory. 4 For discussion of the 
locality of the residual mass term in the transfer-matrix formalism see Refs. lt 

can 

also be seen in perturbation theory 0, 27 



I will use the label m res in various ways in this talk. It has a fundamental, theoretical 



4 The locality of the effective four-dimensional Dirac operator for DWF has been shown to hold (independent 
of L 5 ), under the assumption of sufficiently small lattice field strengths (which forbids zero modes of the 



Wilson kernel), or in the presence of an isolated zero- mode 30(. These are the same conditions under 
which locality of the overlap Dirac operator has been established 3l|. As for overlap fermions, the result 
is not directly applicable to actual DWF simulations, in which no constraint on lattice field strengths is 
applied, and which have a non-zero density of near-zero modes. Thus in practice _one must rely on the 
mobility-edge discussion of Ref . 28j , and check numerically for locality, as in Ref . 
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definition, as a coefficient in the Symanzik effective Lagrangian [eq. (iTTj) below], but with 
this definition it is not directly calculable. It is defined in practice using eq. (|HD below, but 
this includes discretization errors (absent in the fundamental m rcs ) which make it (weakly) 
quark- mass dependent rather than a constant. I follow recent usage and call the calculable 
quantity m f Tes (m) [32, 33[. Finally, I will also use m res generically to mean "a chiral symmetry 
breaking quantity with similar magnitude to the fundamental value" . In this generic usage 
the dependence on L 5 has the same form as in eq. ([2]), but the coefficients of the two terms 
differ from those for the fundamental m res . 5 The main reason for this generic usage is that 
the Symanzik Lagrangian is invalid inside UV divergent loop integrals — such loops give rise 
to quantities with a similar magnitude and form to the fundamental m res , but which are 
numerically different. Furthermore, I will be making order-of-magnitude estimates, involving 
unknown coefficients, so a precise numerical value for m res will not be needed. 

What matters, then, is the approximate magnitude of m res in present and upcoming 
simulations. A lot of work has clarified that if the coupling is too strong (i.e. if a is too 
large) then a becomes too small and chiral symmetry is lost. Happily, one can avoid this 
problem by moving towards weak coupling, and further mitigate it by tuning the action. It is 
also noteworthy that p(0), the coefficient of the power-law term in m res , is expected to drop 
very rapidly as a — > [34]. For present 2 + 1 flavor dynamical simulations at a pa 0.12 fm 
(1/a ~ 1.6 GeV), with L 5 = 16, and the Iwasaki gauge action, one has m res = 0.003 (35|. In 
physical units this is m res /a ~ 5 MeV, i.e. slightly larger than m^ hys . Moving to a ~ 0.09 fm 
(1/a pa 2.2 GeV) at the same L 5 , one finds a further reduction, with m res ^ 0.001 (m rcs /a < 
2 MeV) [36l | . although the precise value is not yet known. 6 Eventually, as a decreases at 
fixed L 5 , one expects from perturbation theory that the dimensionless m res asymptotes to 
a constant up to logarithms, so that m res /a will start to increase. Fortunately, simulations 
have not yet entered this asymptotic regime. 

The bottom line is that in the next couple of years one will have m res pa am^ hys , and the 
question is whether this is small enough to calculate various quantities of interest. 



A. Comparing Wilson and DW fermion analyses of chiral symmetry breaking 



I begin by recalling the classic analysis of XSB for (unimproved) Wilson fermions [37 
since it will serve as the template for the subsequent discussion of DWF. Explicit breaking 
of chiral symmetry leads to a violation of the PCAC relation 

d^Al = 2mP b + aX b , (3) 

where is a particular discretization of the flavor non-singlet axial current, P b the ultra- 
local pseudoscalar density, and m the quark mass (taking degenerate quarks for simplicity). 
aX b arises from the Wilson-term in the fermion action, and thus at tree level has the 
explicit factor of a as shown. This is canceled in matrix elements by factors of 1/a from UV 
divergences, corresponding to the dimension-5 operator X b mixing with lower- dimensional 



5 Note, however, that the coefficient of the exponential, a, does not change (for a given action), since it is 
determined by the mobility edge. I am grateful to Yigal Shamir for emphasizing this to me. 

6 The very recent study of Ref. [8] allows one to determine the relative contribution of the exponential and 
power-law contributions to m' les [see eq. j2j]. One finds that, for £5 = 16, these terms are close to equal 
at 1/a w 1.6 GeV, while at 1/a ~ 2.2 GeV the exponential term is larger by a factor of about 6. 
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operators. To understand the impact of the extra term, then, one must determine the 
operators of dimension four or less with which X b can mix. The answer is that the lattice 
symmetries allow 

777 

aX b ~ -2^P b - (Z A - 1)^4 + 0(a) , (4) 

where ~ indicates equivalence in on-shell matrix elements (so that X b does not come into 
contact with other operators) with p <C 1/a (so that 0(ap) corrections are small). 

The first term on the r.h.s. of (j4j) can be absorbed by an additive renormalization of the 
quark mass, with m^ hjs oc m — m c / a, while the second leads to a renormalization of the axial 
current. After these renormalizations one regains the PCAC relation up to corrections which 
are truly of 0(a), i.e. d^Z^A 1 *) — 2(m — m c /a)P b ~ 0(a). To implement this, however, one 
needs two conditions to determine the two constants m c and One approach is to use 
the vacuum to pion matrix element of eq. (j4j), 

(0|aXV) = -2(m c /a)(0|PV) " ( z a ~ 1)(0|9 m 4|tt) , (5) 

which determines one linear combination of the two constants, 7 while the second condition 
is obtained by enforcing the normalization of Ward identities relating 2-point and 3-point 



functions [37jj. The details are not important here — what matters is that there are two 
constants to determine so one needs two conditions. 

Now consider the corresponding analysis for DWF. The axial current is that generated 
by dividing the fifth dimension into two halves and doing opposite flavor rotations in each, 
and has a lattice divergence [24J 

= 2mJ b + 2J b q . (6) 

Here J5 is the pseudoscalar density built out of fields on the two walls, while J$ q is similar 
in form to J5 but "lives" half-way across the fifth dimension. The J§ q term is due to the 
explicit XSB caused by using finite L5. Thus the situation is formally similar to that with 
Wilson fermions, with J5 corresponding to P (up to differences in normalization), and J§ q 
to aX 8 

In the standard analysis [38[ one now argues that, for low-energy matrix elements, the 
symmetry breaking term can mix only with the pseudoscalar density: 



Jt q -^Jt + 0(a). (7) 



Furthermore, because J§ q resides in the middle of the fifth dimension, this mixing is sup- 
pressed by the decay of the physical modes attached to the walls, so that, unlike the Wilson 



7 The condition is usually written in the equivalent way {0\Z A d ^A^ — 2(m — m c /a)P b \n) = 0, which deter- 
mines Zj[ (m — m c /a). This form makes clear that it is not necessary to determine X b explicitly, and 
that one can implement the procedure for any desired forms of and P b [e.g. with the ultra-local axial 
current rather than the simply local form appearing in eq. {3J]. I write the condition as given in the text, 
eq. (jSJ, as it makes the connection to the subsequent DWF analysis more transparent. 

8 That J$ q does not contain an overall factor of a, unlike aX, is technically correct, but perhaps misleading. 
For finite L5, J§ q does not vanish in the classical continuum limit, whereas aX does. On the other hand, 
as noted above, UV divergences imply that matrix elements of aX do not have an overall factor of a, but 
rather vanish logarithmically (as powers of g(a) 2 ) in the continuum limit. In practice this amounts to 
very little suppression of aX, and is a much smaller effect than the exponential suppression of J§ q . 
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case, m rcs <C 1. The result of this discussion is that one has only an additive renormaliza- 
tion of the quark mass. The single renormalization constant can be determined by a single 
condition, with the standard choice being 38,[3§, 0, 41, 42|: 



a ~ (0|J||vr) ' 1 > 

In other words, one enforces (JTj) as an exact relation in the vacuum to pion matrix element, 
including 0(a) terms. The inclusion of these terms causes m' TeB to depend on m (as we will see 



in more detail below), so it now picks up a prime to indicate that it is not a constant [32|, |33 
It is natural to define a shifted quark mass 



(OI^ANtt) 

m g (m) = m + m res (m) / a = 2 (o\jb\ n ) ' ^ 



which in Wilson-fermion parlance is the "Ward identity mass". Since the r.h.s. is propor- 
tional to m 2 , it follows that the pion mass necessarily vanishes when m q = 0. This occurs 
when m = —m' res (m)/a. 

Comparing the Wilson and DWF analyses, one immediately wonders why the latter is 
more simple, involving one parameter rather than two. Should eq. (JTj) not be replaced by 

J\ q ~ J\ - ^p±d»A\ + 0(a) ? (10) 



The answer is yes — the enumeration of operators is just as for Wilson fermions and eq. fflQ|) 
is the correct form. In practice, however, the standard analysis, and eq. ([7j) in particular, are 
not invalidated. This is because the d^A^ term is highly suppressed, Za — 1 oc m 2 cs < 10~ 5 , 
and can be ignored in practice. 

Why is Z a — 1 suppressed by two powers of m res ? The answer can be deduced from work in 
the literature. Since Zy = 1 for any L 5 , Za — 1 = Za — Zy can be calculated in perturbation 
theory (or using non-perturbative renormalization) by determining the size of the matrix 
elements of the left-handed current (which is associated with one wall) between right-handed 
quarks (attached to the other wall). Thus one needs two crossings of the fifth- dimension, 
one in each direction, to obtain a contribution to Za — 1. (This is the analog of the result 
with Wilson fermions that Za — 1 is proportional to r 2 , with r the coefficient of the Wilson 
term.) In perturbation theory each crossing is exponentially suppressed and proportional 
to m res [26]. Using the transfer matrix argument one reaches the same conclusion in a non- 
perturbative context [29[ — in particular, the structure of the near-zero modes does not allow 
one such mode to cause two simultaneous crossings of the type required for Za — l. 9 

In summary, the fact that some of the terms which enter into a Wilson fermion analysis 
(here Za — 1) receive an extra suppression is an example of how DWF are better than 
'Wilson fermions with an overall suppression of XSB effects" . 



9 The argument of Ref. [29| is made in the context of mixing of left-left four-fermion operators, but it applies 
also to Za — Zy- The observation that the interaction induced by near zero- modes is '"t Hooft-like", and 



thus highly constrained, was made in Ref. 43 1. 
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B. Implications of m rcs ^ in pseudo-Goldstone boson sector 



The analysis so far has absorbed the leading effect of XSB into a shift in the quark mass. 
In this section I recall how this can be seen directly in the Symanzik effective Lagrangian, 
and extend the discussion to sub-leading effects suppressed by powers of a. I consider 
the application to the pseudo-Goldstone boson (PGB) sector, and investigate the relative 
size of errors resulting from XSB (which turn out to be linear in a) compared to (chirally 
unsuppressed) errors of 0(a 2 ). For simplicity of presentation I assume degenerate quarks— 
the generalization to non-degenerate quarks is straightforward. 



The Symanzik effective Lagrangian [44] for DWF has exactly the same form as for Wilson 
fermions j45[, since the symmetries are the same. For on-shell quantities one has 

£s ym . = q(P + Z m m)q + Zm ™ iea q q + ac 5 q(a ■ F)q + 0(a 2 ) , (II) 

where q are quark fields in a regularized continuum theory in which factors of the lattice 
spacing are now explicit. 10 Thus q are proportional (at leading order in a) to the usual 
boundary quark fields of DWF. The factor of Z m is needed since m is the standard lattice 
quark mass [i.e. the same quantity as in eq. (jHj)] and thus must be matched to continuum 
quantities. XSB allows the extra, UV divergent mass term, whose coefficient provides the 
fundamental definition of m res . XSB also allows the Pauli-term, the coefficient of which is 
has the same form (j2j), and is thus also highly suppressed. Note, however, that 

this is only a correspondence of magnitudes, and that the detailed dependence of m res and 
C5 on L 5 will differ, because the ratio of the coefficients of the exponential and power-law 
terms differ. 

We are now in a position to determine the form of the relationship between defined 
in (fill and m' res from (jSJ). To do this we first match £s y m. onto the chiral effective theory, 



following the method of Ref . |47j . This leads to 



C x = ^tr(^S^St) - ^tr(M,S + T$M q ) + ... (12) 



M q = m+ — + ac 5 A 2 , (13) 



This description hides a technical difference between the Wilson and DWF analyses related to the depen- 
dence on a. This point was discussed at the workshop, and here I offer my opinion. The Symanzik action 
for Wilson fermions is of the standard form, with the dependence on a given by the explicit factors of a n 
together with an implicit logarithmic dependence contained in the coefficients of the terms. In the DWF 
case, there is an additional, more complicated, a dependence contained in m ros , eq. ([2]). In particular, as 
noted above, p(0) has a rapid, possibly exponential, dependence on a, and is a non-perturbative quantity. 
Does this invalidate the application of the Symanzik expansion to DWF? I think not, for two reasons — 
one general and one specific to DWF. First, although the Symanzik expansion was originally derived in 
perturbation theory, where one finds powers of a up to logarithms, it can be thought of as an effective field 
theory (EFT), a concept which is not tied to perturbation theory. There is no need for the parameters in 
the EFT to have a simple dependence on a. Second, the presence of the additional parameter L$ allows 
for the appearance of a more complicated a dependence. Indeed, when studying DWF in perturbation 
theory one is effectively constructing terms in the EFT, and one finds a complicated dependence on a, 
with, for example, the exponent a in eq. being a power series in as (a) |27| . 
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where S is the usual non-linear PGB field, / the leading order pion decay constant (nor- 
malized so that f n = 93 MeV, B is proportional to the condensate (and absorbs Z m ) , and 
the dots in (fT2l) indicate terms of higher order in m, m res , a and derivatives. The important 
result contained in eqs. (I12II13P is that the PGB masses come not only from the quark masses 
in £sym. j but also receive a contribution from the Pauli term. This is because the Pauli term 
has the same chiral transformation properties as the quark mass terms. The size of the 
Pauli-term contribution is not, however, known precisely — both because C5 is unknown and 
because the mapping to XPT introduces the scale A which, while of O(Aqcd), is not known 
otherwise. 11 

The net result is that the 0(a) term from Cs ym . is completely absorbed into M q , at 
leading order in XPT. Since m 2 oc M q , and recalling from above that m n vanishes when 
m q = m + m' res /a = 0, one finds 

M q = m q =>• m' rcs /a = m res /a + c^aA 2 , (14) 

up to higher order chiral corrections, and terms of size m 2 es . In other words, when we use 
the standard definition of m' xe8 we are automatically including the contribution to the pion 
masses (and scattering amplitudes, etc.) from the Pauli-term. We do not need to introduce 
an additional XSB term of size c^aA 2 ~ m Tes aA 2 into the leading order C x . 12 

This automatic improvement does not work for higher order terms in the chiral expan- 
sion. The operators qa ■ Fq and qq are different, and their matrix elements not in general 
proportional. The impact of this difference has been worked out to next-to-leading order for 



Wilson fermions [48(, and can be carried over to DWF with appropriate extra suppressions 
of XSB terms. A mnemonic which allows one to see the size of the resulting discretization 
errors is to take the continuum result and make the substitution m q — > m q + m res aA 2 (except 
in the leading order chiral Lagrangian as discussed above). Two examples are 

TTt 2 

— £ ~ U ~ const. [1 + 0{mJA) + 0(m res aA) + 0{a 2 A 2 ) + ...]. (15) 



m q 



Thus XSB gives rise to multiplicative discretization errors of relative size m res aA — i.e. sup- 
pressed by m res compared to the size expected for unimproved Wilson fermions. Similar 
corrections are present in all hadronic quantities. 

Are these errors important in practice, given present values of m res ? I think the answer 
is clearly no. They are numerically tiny and much smaller than the chirally unsuppressed 
0(a 2 ) terms. To make this concrete, consider the situation with present DWF simulations 
at 1/a = 1.6 GeV, and use A = 0.1 - 0.5 GeV. Then 

/0.1 -0.5 GeV\ , , , AS9 

m res aA w 0.003 — — — — — « (2 - 9) X 10" 4 < (aA) 2 w 0.004 - 0.1 . (16) 
\ 1.6 GeV / 

For percent-level accuracy, the XSB-induced errors can be ignored. 



11 While it is redundant in practice to have the product of two unknown quantities, C5 and A 2 , I keep both 
since they originate from different steps in the matching of the lattice theory onto XPT. 

12 As noted in Ref. 47|], the same automatic O(a) improvement of the leading order chiral Lagrangian holds 
for Wilson fermions. The difference here is that the term which is excluded is already suppressed. 
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C. Mass dependence of m' r 



Numerical results for m' res show a weak linear dependence on the bare quark mass m (see, 
e.g., Refs. [H, [33|]). This is understood in the literature to be a discretization error, but 
here I want to determine more precisely its origin and form. Since m' Tes is defined by the 
ratio of the matrix elements of J§ q and J q , one needs to compare the Symanzik expansions 
of these two operators. Combining these, I find (as usual for on-shell matrix elements) 13 



J 5q 



-J b 5 [1 + 0{a 2 m 2 )] + ac b qa ■ F lb T h q + 



(17) 



I have kept operators up to dimension 5, and dropped all contributions proportional to 
m 2 es . This is just eq. flTJ) with the 0(a) term included. Note that m res and C5 are the same 

-this is required in order to reproduce the result ( TH|) for m' res . 

(or, equivalently, linear in m) arises from 
the different mass dependences of (0| Jf|7r) and {0\qa-F r y 5 T b q\ii). The result is a contribution 
quadratic in a, whose relative size is not suppressed by m res : 



quantities as in eq. (TTTT) 



The dominant contribution to m' res linear in m q 



0) [1 + 0(ma 2 A 



(18) 



Another way of stating this is to consider the derivative of \og(m' ms ) with respect to the 
dimensionless bare quark mass: 



1 dm' 



m' ies d(am) 



1 dm' 



-m' les /a 



m' ies d(am) 



0(aA) 



(19) 



rrtasO 



(Note that the derivative can equally well be evaluated at m q 



at the order I 
I find that 



33 



or m = 

am working.) Crudely evaluating this derivative from the plots in Refs. [3 
for 1/a = 1.6 GeV, the l.h.s. is ~ 1.7 for both L 5 = 8 and L 5 = 16. This implies a scale 
A « 3 GeV, which, while very large, is not unheard of. For example, dm 2 /dm q « 5 GeV, 
and this large scale leads to enhanced discretization errors in improvement coefficients for 
improved Wilson fermions (49J. Nevertheless, the large scale does make one worry that 
something is missing from the analysis given here. 

Although the slope is larger than expected, it only produces a small "uncertainty" in the 
value to use for m' Ie8 (m). A natural choice is to use m' Tes (—m' Tes /a), i.e. the value at m q = 0. 
Another commonly used option is m' res (0). The difference between these two is of relative 
size 0(m' Tes aA), a correction of less than a percent even if aA w 2. Because of this, in the 
following I will treat m' Tes as a mass-independent constant. 

A striking feature of numerical results for m' res is that the slope with respect to the valence 
quark mass is significantly larger (by a factor of ~ —2.5) than that with respect to the sea 
quark mass 32, 33|. The corresponding scale, A, approaches 10 GeV. Why such a large scale 
enters is an unresolved puzzle that bears further thought. 



This equation is a hybrid expression in which J| and J| are lattice operators, rather than operators in 
a continuum effective theory, while qa ■ F^T b q is a continuum operator, in which the regularization and 
renormalization is such that no mixing with lower dimension operators is allowed. This equation can be 
obtained by rearranging the Symanzik expansions for the two operators. 
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D. Enhanced m res effects: condensate 

The small size of the generic effects exemplified by eq. (1151) implies that problems 
can only occur if such effects are enhanced. Two enhancement mechanisms are power di- 
vergences, i.e. mixing with operators of lower dimension, and "chiral enhancement", i.e. 
mixing with operators with less suppressed chiral behavior. It turns out that the former is 
much more effective. 

As a first example of enhanced xSB contributions I consider the quark condensate, (qq) . 
Note that one can also calculate the condensate indirectly, using XPT and results for /„., m| 
and quark masses, but it is also of interest to do a direct, XPT independent, calculation. 

To calculate the condensate one needs, in principle, to first take the infinite volume limit 
and then the chiral limit. Both present significant practical challenges, but I assume here 
that these have been overcome. The question I address is whether, after taking these limits, 
and in particular after taking m q — > 0, the result has significant contamination due to the 
lack of exact chiral symmetry. 



This contamination has been considered in Ref. [38|], and I recall and recast the theoretical 
analysis given in that work. Consider the Symanzik expansion of the scalar bilinear, which 
is schematically of the form 

m - m+ x T cs/a + (&) +•••> (20) 

DWF a 2 cont. 

where I am using m res as defined by the Symanzik expansion, eq. fill I) . Here I keep only 
the leading two operators (the identity operator and the condensate itself), and identify 
schematically only the leading terms in the coefficients of these operators. Chiral symmetry 
requires that the coefficient of the identity operator is proportional to the quark mass, and 
thus quadratically divergent, but xSB allows a 1/a 3 divergence proportional to m res . Note 
that since UV momenta dominate the mixing, one cannot use the Symanzik action eq. ffTTj) 
inside loops. Thus m and m Tes /a do not appear in the linear combination which gives m q . 
In other words, instead of x = 1, one expects x = (9(1). 14 Thus, if one extrapolates to 
m q = m + m ms /a = one finds 

lim lim (qq) D WF = (gg)cont + (x - l)—^- + . . . (21) 

m„^0 L^oo a° 



The second term is the result of the finiteness of L 5 . Since x is not known, this term gives 
an uncontrolled error in the condensate. It can be studied and reduced only by increasing 
L 5 — a very expensive proposition. 

How large is the resulting error in the condensate for present simulations? Since (qq) con t ~ 
Aq CD (with a numerical factor not far from unity), the relative error is 

S(qq) m res 



(qq) a 3 A 3 QCD 



(22) 



14 A heuristic discussion of why this arises goes as follows. In the quark loop which gives the condensate, 
higher order terms in the Symanzik expansion (jlip . i.e. the Pauli term, etc., are not suppressed — overall 
factors of a are canceled by factors of momenta k ~ 1/a. These higher order terms are not removed 
by adjusting the bare mass to cancel the contribution from the leading (m res ) term. One reaches the 
same conclusion by comparing expressions for the condensate and for m rcs using the transfer matrix 
formalism [s , 29 1 . 
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This is enhanced by four powers of oAqcd compared to the generic correction exemplified 
by eq. (|T5|) . Numerically, 

5(qq) 3 x 10-3 

<g ff > ~ (0.3/1.6)3 ~ °' 5 (23) 

for the present 0.12 fm simulations, and only slightly smaller for those at 0.09 fm (pa 0.4 if 
m res ~ 10~ 3 ). I conclude that one cannot directly calculate the condensate with present and 
planned DWF simulations, except to check the order of magnitude. 

I stress that these considerations do not apply to the calculations using overlap fermions, 
using which several groups have studied the condensate in small volumes and compared to 
the e— regime predictions of XPT. The chiral symmetry breaking term is absent or negligible 
with overlap fermions. 



E. Enhanced m rcs effects: e'/e 

I now turn to a much more important example: the matrix elements which must be 
calculated to determine whether the Standard Model predicts the measured value of e'/e. 
One needs to calculate the K — > tttt matrix elements of the operators in the AS = 1 
electroweak Hamiltonian, the most important of which is 

Ce = s a 7 M (l - 75K ^2 <?&7m(1 + lb)q a , (24) 

q=u,d,s 

(where a and b are color indices). While ultimately a high precision calculation would be 
desirable, even a result with 10-30% precision would be very interesting. 

As is well known, many challenges must be overcome to calculate these matrix elements. 
One is the difficulty of a direct calculation involving a two pion final state with physical (or 
near physical) kinematics — the method is known in principle, but involves very large lattices 
and the need to pick out excited states [5CJ . This can be circumvented, at the cost of a loss 
of some precision, by using XPT to relate the physical matrix elements to unphysical ones 
which are easier to calculate: 

• At leading order (LO), one can use K — > (m^ 7^ m s ) and K — > tt (m^ = m s ) matrix 
elements [511 ] : 

• At next-to-leading order (NLO), one must add further matrix elements to determine 
all the needed low-energy constants. One scheme uses K — > tt with 7^ m s , K —>■ K 
matrix elements of the AS = 2 Hamiltonian, and unphysical K — >• tttt matrix elements 
with the pions at rest [52{|. 15 Another uses K — >• tttt matrix elements involving particles 



with non-vanishing momenta 53 . 



The precision that one can obtain using these two methods is unclear. Naively, for kaon 
matrix elements, the relative errors in the LO calculation are of size (mjf/A x ) 2 pa 0.2, where 
A x = 47t/, while those in the NLO calculation are ~ (m^/A x ) 4 pa 0.04. There is considerable 
evidence, however, that these errors are larger for e', and it may require a NLO calculation 



15 



One can also use partiall y q uenched matrix elements to provide further information and cross-checks. 



This is discussed in Refs. 



43 



52| . For simplicity, I do not discuss this extension here. 
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to achieve the desired 10-30% precision. In practice, the LO calculation will be done first 
(and is indeed underway now), with the NLO refinements added subsequently (5^]. The LO 



methodology has been worked out and tested in quenched calculations [46j, [55 



This brings me to the point I want to address here — how is this calculation affected by the 



lack of exact chiral symmetry? This has been discussed in part by Refs. [46j, [55j, and more 
extensively in Ref. |43j ]. but here, based on considerations with Wilson-like fermions (BH 57], 
I extend the discussion. 

Even with exact chiral symmetry, as soon as one considers unphysical matrix elements 
one must account for mixing with the lower-dimension operator |5l| 16 

Csubi = \{(m d - m s )sj B d + (m, + m d )sd} . (25) 

In the XPT-based methods discussed above this is subtracted using the K — » matrix 
element. This only works, however, if the positive and negative parity parts of the operator 
are related as in O su bi> which requires chiral symmetry. The breaking of chiral symmetry 
with DWF leads to additional operators to subtract. I list here those of dimension less than 
6 (which thus have power divergent coefficients and give enhanced contributions) which are 
linear in m res (as opposed to quadratic or higher order) 

Osub2 = ^sd, (26) 
Osub3 = ^StT-Fd, (27) 

0£m = —M, (28) 

m res (mg - m 2 d ) 

(J sub 5 = S75GL (29) 

a 

Because of the presence of UV divergences, it is the generic m rPK that appears in these 



operators. This list is almost the same as for Wilson fermions 56], except that here each 



occurrence of a 11 L <-> R flip" is always suppressed, either by a factor of m res or of m q . Note 
that the operators Of nhA and O su b5 are usually subsumed into mass-dependent coefficients of 
O su bi and C su b2- I prefer to keep the new operators explicit, since the factors of m q impact 
the fits to XPT. 

The presence of these operators invalidates the XPT-based methods, and if not sub- 
tracted, they lead to errors in matrix elements that can only be reduced by increasing L 5 . 
Clearly the greater the power of the UV divergence, the larger the potential error. 

The presence of the most divergent operator, O su b2, was noted in the quenched stud- 
ies (46, 55|, and a method was proposed and implemented to subtract it within the LO XPT 



approach. I recall this method briefly here as I will build on it below. In order to determine 
(K\Oq\tttt) at LO in XPT one needs to calculate, for m s = ma, 

M phys = (K\O e - ciO subl - c 2 O sub 2|vr) , (30) 



16 Here and in the following I will assume for definiteness that the quark masses include the shift due to mj esl 
i.e. are given by eq. ([9]), although one could equally well use bare lattice quark masses for this analysis. 

17 In fact, this is not quite correct: near-zero modes allow multiple flips at the cost of a single m res sup- 
pression [43j |. As explained in the talk by Christ 58|, this allows an operator of the form (m res / a 2 )m q sd, 



which would naively require an additional factor or m ros . I discuss this operator further below. 
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where I have dropped O su b3 etc. for now. The subtraction coefficient c\ is to be determined 
from the K — > matrix element, and I assume that this can be done with sufficient accuracy, 
as seems likely from the quenched studies. The problem then arises because the coefficient C2 
is not known. Without a determination of C2, the matrix element will have an unacceptably 
large error 



5M phys C2 TTLres 

B/a 3 Aqcd 



(31) 



•Mphys m 2 K A 2 QCB (aA QCD ) 3 m s 

Here I use the estimate -M p h y s ~ m |-^-QCD; which incorporates the chiral suppression and 
is consistent with numerical results from quenched simulations and continuum estimates. 
I also use the LO XPT result {K\sd\n) « B « m 2 K /m s , and take c 2 ~ 1, m s « 90MeV, 
Aqcd = 300 MeV, and, finally, m res = 0.003 for 1/a = 1.6 GeV, as above. Note that the 
errors are of the same form as the condensate, eq. (1221) . except for the additional chiral 
enhancement. 

This error due to C2 7^ can be removed, however, with a second subtraction condition. 
One approach is to enforce that -M p h y s vanishes in the chiral limit, since this holds in a 
chirally symmetric theory. The idea is then to adjust C2 until this vanishing holds. If 
one works at LO in XPT it is then the slope versus m s = that gives the required 
information [46|, 55[. At NLO a more sophisticated fit to the data is required, including 



quadratic terms and m s 7^ m<i, but the same subtraction condition works. 18 

If O su b2 can be subtracted, then the largest uncertainty is likely due to the next-most 
divergent operator O su ^. The presence of this operator is well known to Wilson and DWF 
practitioners, but it is only now, with the advent of realistic DWF simulations, that it 
becomes relevant to estimate its contribution. Since its K — > it matrix element does not 
vanish in the chiral limit, it contaminates the subtraction method for (9 su b2 — a single con- 
dition (vanishing in the chiral limit) cannot determine both unknowns c 2 and C3. One must 
either accept the resulting error or determine another subtraction method. The former ap- 
proach may, in fact, be sufficient. One expects that the ratio of matrix elements of sa ■ Fd 
and sd is approximately Aq CD . Using this, and the parameters given above, one finds 



5-Mphys c 3 m res A QCD m r 



-Mp hys aA QCD m s am s 



0.05 (33) 



for the a = 0.12 fm simulations. At a = 0.09 fm, the result is 2-3 times smaller. Thus the 
error one makes by ignoring this term is comparable to, or smaller than, that due to the 
missing next-to-next-to-leading (NNLO) contributions in the XPT analysis. It seems to me 
that one can, at this stage, proceed by ignoring this operator. 

It should also be possible to subtract the contributions of this operator. One idea (which 
I have not thought through in detail) is to use the fact that -M p h y s is proportional to m^rriK, 



One can avoid the need to rely on accurate fitting of the matrix element by using the ratio 

(K\Q 6 - ci0 sub i|7r) _ A^phys ra rcs , . 

(K\sd\ir) ~ (K\sd\w) a 3 ' ( ' 

The first term on the r.h.s. has a chiral expansion beginning at linear order in m q , so it should be 
straightforward to determine and remove the C2 term, which is a constant. One can then reconstruct 
Afphys by multiplying by the separately measured scalar matrix element (K\sd\ir). This method is due 
to the RBC collaboration. 
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whereas the matrix element of (9 S ub3 also has terms of the form m 2 and m 2 K (which are non- 
leading for this operator but will contribute to the putative -M p h ys if C S ub3 is not subtracted). 
The idea would then be to adjust C3/C2 so as to remove these terms. Another approach is 
suggested in Ref. 56|. Both ideas are valid only at LO in XPT and would need to be 



elaborated in order to work at NLO. 

I now return to the remaining operators in the list, O suhA ± and S ub5- Given the explicit 
factors of quark mass, these do not affect the LO XPT analysis, but do contaminate the 
NLO analysis. They lead to an error which is suppressed by (am s ) 2 relative to that which 
would be present if (9 su b2 were not subtracted [see eq. fl31~]) ]. and thus of size 



5M ph y S C 4i5 m res Aq C D^_ a2 m res ™? 

V QCD 



— — - K am s ) 2 ~ — -ry^- ~ 0.005 . (34) 

.Mphys (aAQCDr m s am s A( 



It appears, therefore, that one can safely ignore these operators for a long time to come. 

The discussion thus far turns out to be incomplete — there is an additional contribution 
noted by Norman Christ in his talk here j58[ . This arises from an operator which is naively 
be suppressed by m 2 es m q , but is in fact suppressed only by m res m q because a single near- 
zero mode can simultaneously flip the chirality of two quarks as long as they have different 
flavors (and the chirality flips are in the "same direction"). The fact that the zero-modes 
are localized causes no suppression because the initial operator itself is local. 19 Christ shows 
that one obtains the operator 

{{m s - m d )sj 5 d + (m s + m d + 2m u )sd} . (35) 

Part of this operator is removed by the subtraction for (9 su bi, allowing one to rewrite it as 

su b6 = ^^2(m s + m d + m u )sd . (36) 
cr 

How does the error due to this operator compare with that from O su b3? There are 
competing factors here. On the one hand C S ub6 is enhanced by one power of 1/a, but, 
on the other hand, it is chirally suppressed by one factor of m s . Naively combining these 
factors one obtains an enhancement of m s /(aAQ CD ) « 2. There is another suppression 
factor, however, namely that it is only the near-zero mode part of m res , i.e. the second term 
in eq. (J2J), which contributes. As noted above (footnote 5) the near-zero mode contribution 
to m! res is, for present simulations, at most comparable to that from the extended modes. 
Thus it seems plausible that the m res appearing in (|36|) is smaller by at least by a factor of 
2 than that appearing in the error estimate due to C su b3. All in all, I expect that C su b3 and 
C su b6 give comparable errors, and thus, based on the estimates given above, that both can 
probably be ignored at this stage of the numerical calculations. 



Eventually, however, one will want to subtract O su b6 as well. As noted in Ref. [58|], this 
is possible, in principle, by varying m u independently from m d and m s . 

In conclusion, as long as one implements a method to subtract the most divergent operator 
mixing caused by xSB, DWF can probably be used to study e'/e using XPT-based methods 
(particularly the NLO method) with sufficient accuracy to make significant progress. 



19 This is in distinction to the generic contributions from the near-zero modes, which, as explained in Ref. 43 1 
are suppressed by additional factors of a due to the localized nature of the eigenmodes. 
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F. Implications of m rcs for Bk 

I now turn to another phenomeno logically important matrix element, Bk- This is an 
example where there is no mixing with lower dimension operators, but mixing with chirally 
enhanced operators can occur. Thus one might expect a "chiral enhancement" of XSB 
effects. It is also a quantity where percent level precision is highly desirable, since the better 
we can calculate Bk-, the more precisely we can test the Standard Model. 

Chiral enhancement can occur through mixing of the original Bk operator, Ok = 
SLlfidLSL'JfidL, with "L-R" operators exemplified by SLj^dLSjij^dR. The L-L operators have 
K — K matrix elements which are chirally suppressed, and thus proportional to m s , while 
the L-R operators are unsuppressed. Since two quarks must have their chirality flipped, 
mixing at the perturbative level is proportional to m 2 es (the argument being identical to 
that for Za — 1 discussed above). 20 Even with a chiral enhancement of size Aqcd/^s, this 
contribution is numerically insignificant. Including discretization errors [i.e. working out 
the Symanzik expansion for Ok arid keeping terms analogous to the Pauli-term in (111}) ] I 
find that there is a mixing proportional to m res am s . The factor of m s is needed because of 
chirality — one flip requires am res and the other m s . Here m res is meant generically, since the 
contribution involves logarithmically divergent integrals and thus UV momenta. The extra 
factor of m s in the mixing cancels the chiral enhancement in the matrix element, so one 
finds in the end an error which has the same size as those in generic hadronic quantities, i.e. 

Aqcd . /o n x 1A _4 /0 „n 

— — ~ m ves am s — — ~ m res aAQCD ~ (2 - 9) x 10 . (37) 
B K m s 

Here the numerical estimate is from eq. ({TBI) . Errors of this same size also arise directly from 
the Pauli-term in the Symanzik Lagrangian. 

I conclude that XSB is not a concern for percent level calculations of Bk with DWF. 
Indeed, such precision is already being attained at a single lattice spacing [4J , with a second 
lattice spacing underway (H3 |. 



G. Pion EM splitting 

I close with a final example which is both of phenomenological interest, and has particular 
resonance in this meeting as it was based on a suggestion (long ago) by David Kaplan. The 
idea is to calculate the pion EM mass splitting using XPT to relate it to a simpler quantity: 

2 2 /" 1 

A. = =-J ^i^^AOr) d Rl ,u R (0)) , (38) 

where the r.h.s. should be extrapolated to the chiral limit. The r.h.s. is part of the vacuum 
energy due to electromagnetism, with the 1/x 2 arising from the photon propagator. This 
application of (what was then called) current algebra and PCAC was suggested in 1967 by 
Ref . [59] , and they evaluated the correlation function using resonance saturation constrained 
by the Weinberg sum rules. For an updated discussion incorporating large N c input, see, 



e.g., Ref. [60 



21] 



Note that near-zero modes do not give a contribution proportional to m res here, just as for Za — 1 2S|. 
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Of course one can also do a direct calculation of EM splittings, using a background EM 
field, and indeed such calculations are underway using DWF j6l|. But it would be interesting 
to have a result using a different method, and also to study correlation functions like this 
so as to make contact with large- N c based approaches discussed here by Golterman (62| . 

A first lattice calculation (which happens also to be my first lattice calculation) was 
attempted long ago, using Wilson fermions on a lattice of size too small to mention [63]. 
We understood at the time that chiral symmetry was essential for evaluating the r.h.s. of 
fl38|) — it is a chiral order parameter, like the condensate, since its L-R structure means that 
it vanishes in perturbation theory to all orders for massless quarks. It quickly became clear 
that the inherent to Wilson fermions introduced an intolerable error. Given the much 
improved chiral symmetry of DWF, I thought it interesting to revisit this quantity and see 
how well one can do. 

To study this, I note that, if m res = 0, one expects, schematically, that 

K~ ft + ^t + m q {m)^a + ... , (39) 
cr 

where the first term is the desired result coming from long-distance contributions to the 
integral over x (with being a good estimate of its magnitude using the experimental 
pion splitting), while subsequent terms are UV divergent and can be determined using 
the operator product expansion. (In doing so, I have not kept track of factors of 0:5 in 
coefficient factors, since I am making rough estimates.) As in the case of the condensate, 
one can determine the order of magnitude of the contribution due to m res by substituting 
m q — > m res /a. In this way one finds that, after chiral extrapolation, the relative error is 

^-f^ + ^f X<§>~ (1) + 0(1). (40) 

In words, both contributions quadratic and linear in m res turn out, for present and near- 
future DWF simulations, to give errors of order 100%. I conclude that it is probably not 
practical to calculate in this way, but it is perhaps worth a more detailed look. 



IV. CONCLUSIONS 

It appears that the residual breaking of chiral symmetry is problematic only for highly 
UV divergent quantities: in the examples I have considered, these are the condensate and 
the pion EM mass splitting. Thus I take the middle path and conclude that DWF are 
WONDERFUL: there is no practical barrier to calculating many quantities of interest 
(spectrum, Bk and related matrix elements, nucleon properties, and possibly e'/e) 
with desired precision in the next five years. This will require, however, simulations with 
parameters extending to at least me/m s ~ 0.1 and L « 4fm, and probably down to a ~ 
0.06 fm, while keeping m res « 10 -3 . Given present estimates of CPU requirements, these 
appear to be attainable parameters in a five year time-scale. 

Given the exciting potential of DWF simulations, I cannot resist an exhortation. For an 
outsider, a striking lacuna in present DWF studies is the lack of calculations of quantities 
involving heavy (and in particular b) quarks. There are so many interesting and important 
heavy-light quantities to calculate that this is a serious omission. More generally, as stressed 
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here by Soni, the lattice community now has a great new method, and needs to think hard 
about increasing its repertoire. 

Clearly, the next such DWF meeting should be very interesting! 
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